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ABSTRACT 

A  synthetic  sampling  frequency  has  been  developed  utilizing  a  25 
kHz  SERVO-channel  CW  signal  and  a  predetermined  frequency  of  choice  to 
achieve  the  proper  sampling  frequency  for  quadrature  demodulation. 

The  required  synthetic  sampling  frequency  examined  here  is  Af  n  and  is 
compared  with  the  actual  pre-recorded 
shown  for  signals  which  were  envelope 
demodulation  methods.  This  technique 
sampling  frequency  of  choice,  and  can 
signal  processing  tool. 


4f  .  Very  favorable  results  are 
detected  utilizing  quadrature 
may  be  applied  to  produce  any 
thereby  prove  to  be  a  valuable 
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INTRODUCTION 


In  the  process  of  A/D  conversion,  one  is  always  faced  with  the 
task  of  choosing  the  appropriate  sampling  frequency  (f  )  so  that  the 
digital  form  of  the  data  is  unbiased  and  truly  reflects  its  analog 
counterpart.  This  is  all  so  apparent  when  applying  the  techniques  of 
quadrature  demodulation  to  the  envelope  detection  of  a  signal.  With 
one  technique  of  quadrature  demodul at  ion,  consecutive  samples  taken 
exactly  90°  apart  at  a  sampling  frequency  of  4f  (four  times  the 
center  frequency),  are  needed  for  the  demodulation  to  be  accomplished 
correctly.  A  more  detailed  explanation  of  this  quadrature 
demodulation  technique  is  given  in  Appendix  I. 

During  the  A/D  conversion  process  of  some  recently  acquired  analog 
tapes  from  the  NOREX  '83  experiment,  a  problem  arose  when  the  tape 
recorded  version  of  the  4f  was  missing  from  some  of  the  analog  tapes. 
The  engineering  task  was  t8  develop  a  synthetic  4f  which  had  the  same 
attributes  as  the  real  4f  and  could  be  used  to  su2cessfully 
quadrature  demodulate  the°data  signals.  As  will  be  shown,  there  is  a 
relatively  straightforward  solution  to  the  task  as  long  as  a  CW  signal 
has  been  recorded  which  preserves  the  phase  errors  introduced  by  tape 
speed  variation. 

The  hardware  set-up,  and  subsequent  testing  of  the  synthetic  4f 
against  a  real  4f  is  presented  in  the  sections  which  follow. 

HARDWARE  SET-UP 


The  hardware  set-up  for  the  generation  of  the  synthetic  4f  is 
shown  in  Figure  1.  For  this  case,  the  25  kHz  SERVO-frequency  on 
analog  tape  is  band-pass  filtered  and  fed  into  a  multiplier  box.  A 
predetermined  mixing  frequency  of  constant  phase  is  generated  by  a 
frequency  synthesizer  and  inserted  into  the  mixer  [see  Appendix  II, 
case  II].  The  output  of  the  multiplier  consists  of  two  frequencies,  a 
sum  and  a  difference  frequency.  The  output  of  the  multiplier  is  then 
bandpass  filtered  for  the  desired  sum  or  difference  frequency. 


For  example,  if  one  is  interested  in  a  quadrature  demodulated 
signal  with  a  center  frequency  (f  )  of  20  kHz,  a  sampling  frequency 
(f  )  of  4f  ,  or  80  kHz  is  needed  in  real  time.  If  this  is  done  at 
l/8th  real°time,  as  was  the  case  with  the  NOREX  '83  20  kHz  data,  a  10 
kHz  sampling  frequency  is  actually  required.  In  this  case,  the  25  kHz 
SERVO-frequency  from  tape  is  at  3.125  kHz  and  the  predetermined  mixing 
frequency  is  6.875  kHz.  The  result  is  a  sum  frequency  of  10.0  kHz  and 
a  difference  frequency  of  3.75  kHz.  This  complex  signal  is  then 
bandpass  filtered  at  10  kHz  to  achieve  the  desired  4 f  . 


The  synthetic  4f  was  then  compared  with  a  real  4f  Q  ,  which  was 
recorded  on  tape,  to°determine  the  validity  of  the  approach,  as 
illustrated  in  Figure  2. 
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RESULTS 

A  comparison  of  the  real  4f  and  the  synthetic  4f  was  done  for  a 
signal  with  a  center  frequency  of  20  kHz.  For  the  case  examined,  the 
same  processing,  quadrature  demodulation,  and  A/D  conversion 
procedures  were  used  on  the  same  data  set  of  CW  pulses,  the  only 

difference  being  the  type  of  4f  sampling  frequency  used  in  the  A/D 

conversion  process.  In  one  case  the  real  4f  from  tape  was  used  [see 

Appendix  II,  case  I],  and  in  the  second  case  a  synthetic  4f  was  used 
as  described  in  Appendix  II,  case  II.  The  comparison  is  made  on  an 
individual  pulse  basis  and  also  on  the  ensemble  average  of  twenty-five 
pulses. 

In  Figures  3a  and  3b,  the  level  versus  time  of  the  signal  envelope 
is  presented  for  comparative  purposes,  for  the  real  and  synthetic  4fQ 
sampling  frequencies.  It  is  readily  apparent  that  on  a  pulse  for 
pulse  and  ensemble  average  basis,  there  is  no  observable  difference  in 
the  pulse  or  signal  envelopes.  Due  to  A/D  noise,  differences  are  seen 
in  noise  values  as  expected  [see  Appendix  III], 

Since  the  area  of  interest  is  the  signal  envelope,  Figure  4 
represents  an  enlarged  version  of  a  representative  envelope  which  is 
presented.  Here,  a  well  defined  smooth  envelope  is  seen  for  the  real 
4f  and  synthetic  4f  .  If  the  quadrature  demodulation  technique  was 
no?  being  accomplished  properly,  a  decrease  in  level  or  a  modulation 
component  would  be  seen  on  top  of  the  envelope  [see  Figures  A1  and 
A3].  Since  neither  is  occurring,  we  can  conclude  the  quadrature 
demodulation  technique  is  being  implemented  properly  with  a 
synthetically  generated  4fQ. 

With  the  magnitude  portions  of  the  real  and  synthetic  4f  in 
excellent  agreement,  the  final  aspect  to  make  the  comparisons  complete 
is  the  phase  of  the  signal.  Figures  5a  and  5b  present  a  pulse  to 
pulse  comparison  between  the  two  types  of  sampling  frequencies  for 
four  different  pulses.  Once  again  slight  differences  are  seen  in  the 
phase  of  the  noise  floor,  but  the  phase  components  of  the  signal  are 
in  excellent  agreement. 

With  excellent  agreement  in  magnitude  and  phase,  a  check  was  done 
on  an  ensemble  of  pulses.  Figure  6  represents  this  comparison  and 
once  again  excellent  agreement  is  seen  between  the  real  4f  and 
synthetic  4fQ  with  only  0.1  dB  difference  in  the  mean  values. 

With  the  completion  of  this  last  comparison  a  summary  is  now 
presented  in  Table  1.  Here,  very  good  agreement  is  seen  in  the  first 
order  statistics  between  the  real  and  synthetic  4f  sampling 
frequencies.  The  largest  differences  in  this  table  are  once  again 
associated  with  the  noise  floor. 
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SUMMARY  AND  CONCLUSIONS 

A  4f  sampling  frequency  was  required  in  order  to  implement  a 
quadrature  demodulation  technique.  There  are  however,  no  restrictions 
which  preclude  extending  the  results  of  this  work  to  other  sampling 
f requenci es . 

After  some  theoretical  considerations  and  experimentation,  a 
synthetic  4f  was  produced  which  did  possess  phase  and  frequency 
characteristics  similar  to  the  real  4f  .  This  synthetic  4f  was 
produced  by  multiplying  a  25  kHz  CW  tape  recorded  SERVO  signal  and  a 
predetermined  mixing  frequency.  The  25  kHz  CW  SERVO  signal  from  tape 
had  the  same  time  dependent  phase  as  the  signal  of  interest,  and  the 
predetermined  frequency  CW  from  the  frequency  synthesizer  had  a  time 
independent  phase.  The  multiplier  output  was  subsequently  filtered  to 
discriminate  between  the  sum  and  difference  sampling  frequency. 

A  set  of  data  were  processed  using  a  recorded  (real  4f  )  and 
synthesized  (synthetic  4f  )  sampling  frequency  to  form  comparative 
test  cases.  The  results  8 f  the  comparison  between  amplitude,  phase, 
arrival  time,  and  first  order  statistics  of  the  signal  envelope  were 
excellent. 

Since  the  results  between  the  synthetic  4f  and  the  real  4f  are 
in  such  good  agreement,  it  is  my  opinion  the  synthetic  4f  process  can 
indeed  be  used  for  the  A/D  conversion  process. 

If  this  synthetic  sampling  frequency  technique  is  used, 
considerable  thought  should  be  given  to  selecting  the  CW  SERVO  signal 
which  will  be  recorded.  This  is  a  necessity  because  when  the  CW  SERVO 
signal  is  subsequently  multiplied  and  filtered  to  obtain  the  synthetic 
4f  ,  the  sum  and  difference  frequency  are  far  enough  apart  so  that  a 
filtering  problem  is  not  encountered  in  trying  to  segregate  the  sum 
and  difference  frequency. 

If  this  technique  is  to  be  implemented  without  unnecessary 
problems,  possible  drawbacks  of  the  technique  should  be  highlighted. 
The  first,  and  most  important  drawback  is  the  choosing  of  a  proper  CW 
SERVO  signal  on  tape  so  as  not  to  impose  filtering  problems  when  this 
signal  is  multiplied  with  the  mixing  frequency.  Another  drawback  is 
that  the  actual  production  of  a  synthetic  sampling  frequency  will 
entail  utilizing  more  equipment  than  required  if  the  desired  sampling 
frequency  were  recorded. 

In  closing  it  should  be  noted  that  the  author  highly  recommends 
this  technique  of  producing  a  synthetic  sampling  frequency  for  use 
when  the  optimum  sampling  frequency  has  not  been  pre-recorded. 
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APPENDIX  I 

OVERVIEW  OF  A  QUADRATURE  DEMODULATION  TECHNIQUE 

Let  x.(t)  be  a  signal  from  hydrophone  i  and  y.(t)  be  the  90°  phase 
shifted  signal  from  hydrophone  i. 

The  time  sampled  signals  can  be  represented  as;  x.  (mAt),  y.  (mAt) 
where  m  is  the  digital  sample  point  and  At  is  the  inverse  of  them 
sampling  frequency. 

1,  L 

The  envelope  of  the  signal  from  the  iL  hydrophone  can  be  obtained 
utilizing  the  following  equation; 


As  an  example,  let 


x . ( t )  =  Bsin(wt) 


(1) 


(2) 


y.(t)  =  Bsin(wt  +  90°)  =  Bcos(wt) 


(3) 


where  B  is  constant. 

The  envelope  of  the  signal  [|S.(t)|]  can  then  be  obtained  as 
f ol 1 ows : 

|S,(t)|  =  Cx1 2(t)  +  y,2(t)]1/2  (4) 

|  S ,  ( t )  |  *  [B2($in2wt  +  cos2wt)]^2  =■  B  (5) 

Similarly,  the  phase  of  the  envelope  may  be  defined  as: 


♦1  =  tan  "  1[xi/yi]  (6) 

We  see  there  is  no  reduction  in  the  amplitude  of  the  original 
signal  when  the  signal  is  envelope  detected  using  this  quadrature 
demodulation  technique. 

We  can  achieve  this  result  if  the  sampling  frequency  chosen  fits 
the  criteria  below: 

(center  frequency  (f  ))  (360°)  =  90°  (7) 

(sampling  frequency  ( f  g  ) ) 

For  this  to  be  true  the  sampling  frequency  f  must  equal  4f  (four 
times  the  center  frequency).  s 
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Variations  from  4f  will  cause  an  error  in  the  true  amplitude 
level  of  the  envelope  8s  shown  in  the  example  below. 

Let , 

f  =  20  kHz 
f°  =  81  kHz 

then  from  Eq.  7;  (20.0  kHz)  (360°)  =  88.9 

(81.0  kHz) 


A  88.9°  phase  shift  is  seen  thus  the  identify  transformation  of  (3) 
cannot  be  employed  and  the  term  in  parenthesis  in  (5)  is  not  identical 
to  1.  Therefore,  an  amplitude  error  in  the  envelope  is  introduced 
when  samples  not  spaced  90°  apart  are  used. 
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APPENDIX  II 

Phase  Relations  Between  Signal  and  Sampling  Frequency 
f  -  center  frequency  of  interest 

x. (t)  -  signal  from  a  hydrophone  i 

y .  ( t )  -  signal  delayed  by  90°  from  a  hydrophone  i 
f  -  sampling  frequency  =  4f 

<(>  -  phase 

CASE  I  -  Hydrophone  signal  and  correct  sampling  frequency  on  tape. 


xi(t)=Asin[w1t+«fr1(t)]  (1) 

f  s  ( t )  =  Bsi  n[4w^t  +  «^(t)]  (2) 

where  the  quantity  in  brackets  [  ]  will  be  called  0.  The 

instantaneous  value  of  0  for  the  above  two  expressions  is 

Q[x  i  ( t )  ]  =  Wj  +  <{>i  (  3  ) 

e[  f  s  ( t )  ]  =  4Wj  +  ^  (4) 


where  the  dot  represents  the  time  derivative.  It  is  apparent  the 
phase  terms  tof  the  signals  exhibit  the  same  time  dependence,  <|>*  from 
(3)  equals  <|>*  from  (4),  therefore  quadrature  demodulation  can  oe 
accompl ished1successfully. 

CASE  II  -  Signal  and  CW  SERVO-f requency  on  tape. 

The  sampling  frequency  is  achieved  by  multiplying  the  CW  SERV0- 
frequency  by  a  properly  predetermined  frequency  from  a  frequency 
synthesizer  the  phase  of  which  is  time  independent.  The  phase  term 
from  any  signal  on  tape  is  time  dependent  because  of  tape  speed 
variation. 

Here,  the  multiplier  output  is 
[si  n(  w  t  +  <(>  (t)]sin(w  t  +  <k) 

i  2  2  1 1  3  .  i 

i  i 

CW  SERVO  frequency  synthesizer  frequency 
and  the  result  obtained  is 

l/2[cos(w2t  +  <|>2  ( t )  -  w3t  -  <{>3)  -  cos(w2t  +  <|>2(t)  +  w3t  +  4»3 )  ] 

=  1/2[cos((w2  -  w  3 )  t  +  <f2(t)  -  <f3)  -  c  o  s  (  ( w  2  +w3)t  +  <t>2(t)  +  4>3 )  3 
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After  filtering  the  undesired  frequency  we  obtain  the  proper  4f  .  The 
signal  and  sampling  frequency  are  now  shown  below. 

xi (t)  =  Aslntwjt  +  ♦j(t)]  (5) 

fs(t)  =  Ccos[(w2  -  w 3 ) t  +  <t»2(t)  -  <t>3]  (6) 

As  in  case  I,  our  concern  is  in  the  instantaneous  values  of  the  above 
functions  of  9  (terms  in  brackets)  for  equations  (5)  and  (6); 

9[x1  (t )]  =  Wj  +  (7  ) 

9[f  5  ( t )]  =  w2  -  w3  +  {2  =  w2  -  w3  +  {j  (8) 

where  w2  -  w3  =  4w,.  Therefore,  equations  (4)  and  (8)  are  equivalent 
and  once  again  quadrature  demodulation  can  be  used  successfully 
because  the  signal  and  instantaneous  sampling  frequency  exhibit  the 
same  time  dependence. 


CASE  III  -  Only  the  signal  is  on  tape,  and  sampling  frequency  from  a 
frequency  synthesizer. 


x. (t)  =  As1n[Wjt  +  ♦j(t)] 
f$1(t)  =  Dsi  n[4wxt  +  <t>3] 
f  s  2  ( t )  =  Dsi  n[4Wjt  +  <t>3(t)] 
Taking  the  instantaneous  value  of  9 
9(x. (t ) )  =  w l  + 

9(f$1(t))  =  4wj 
0(fs2(t ) )  =  4Wj  +  i3 


(9) 

constant  phase 

(10a) 

varying  phase 

(10b) 

again  reveal s ; 

(ID 

(12a) 

(12b) 

It  is  apparent  there  is  a  phase  difference  between  (11)  and  (12a),  and 
also  between  (11)  and  (12b).  With  this  difference  in  time  dependent 
phase  between  the  signal  and  the  sampling  frequency,  quadrature 
demodulation  cannot  be  implemented  successfully.  The  only  way  this 
could  occur  is  if  <L(t)  =  (t).  This  is  difficult  to  perceive  since 

the  signal  and  the  Sampl  i ng frequency  are  from  two  physically 
different  pieces  of  equipment. 


In  Figures  A1  through  A4 ,  the  importance  of  the  same  time 
dependent  phase  between  signal  and  sampling  frequency  is  illustrated. 
In  Figures  A1  and  A2 ,  the  differences  are  very  evident  on  a  pulse  for 
pulse  basis.  We  see  a  difference  in  level  as  well  as  modulation  on 
top  of  the  envelope  of  the  signal  in  Figure  Al.  In  Figures  A3  and  A4 
we  can  see  the  differences  in  the  ensemble  averaged  envelope  that  is 
not  apparent  in  Figures  Aid  and  A2d.  Here,  the  presence  of  a 
modulation  on  top  of  the  envelope  is  also  seen. 
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APPENDIX  III 
EFFECTS  OF  A/D  NOISE 


In  situations  of  high  signal-to-noise  ratio,  the  lower  limit  of 
the  A/D  converter  is  rapidly  reached  and  the  resultant  values  may  be 
subject  to  noise  error.  In  this  case,  a  12  bit  A/D  was  used  with  a 
dynamic  range  of  66dB.  This  means  that  if  a  data  set  with  a  high 
signal  to  noise  ratio  is  digitized  two  independent  times,  one  expects 
to  see  exact  replication  of  the  signal  envelope,  and  some  random  error 
introduced  into  the  low  noise  values  due  to  the  A/D  process. 

In  Figure  A5  the  same  data  set  is  digitized  twice  with  the  real 
4  f  .  The  exact  duplication  of  the  signal  envelope  is  seen,  along  with 
the  non-exact  replication  of  the  noise  floor.  Figures  A6  and  A7 
illustrate  the  noise  level  of  the  A/D  system  which  is  shorted  at  the 
data  input  and  shorted  at  the  input  to  the  computer,  respectively  (see 
Figure  2).  It  becomes  apparent  that  the  noise  levels  due  to  the  A/D 
system  (filters,  amplifiers,  and  computer)  are  nearly  the  same  as  the 
noise  levels  of  the  computer  A/D. 

Therefore,  if  noise  floor  level  random  errors  are  introduced  when 
the  real  4f  is  utilized,  they  would  surely  be  expected  when  the 
synthetic  4r  is  used. 
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Figure  1.  Basic  hardware  set-up  for  the  production  of  a  synthetic  sampling  frequency. 
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Figure  2.  Production  of  a  synthetic  4f0  vs.  the  real  4fo  for  a  20  kHz 
center  frequency  (system  is  at  1/8  real  time). 
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Figure  4.  Enlargement  of  envelope  of  one  pulse  utilizing  a)  real  4fo  and 
b)  synthetic  4fo. 
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lses  utilizing  the  real  4fQ  sampling  frequency. 
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(a)  (b) 

Figure  6.  Peak  envelope  level  of  an  ensemble  of  25  pulses  utilizing  (a)  real  4fo  and 
(b)  synthetic  4fo>  for  the  quadrature  demodulated  envelope. 
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CRITERIA 

REAL 

4*o 

SYNTHETIC 

NOISE 

(dB) 

MEAN 

-32.42 

-32.78 

STANDARD 

DEVIATION 

1.24 

1.45 

COEFFICIENT 
OF  VARIATION 

0.038 

0.044 

LEVEL 

OF 

SIGNAL 

ENVELOPE 

(VOLTS) 

MEAN 

3.38 

3.32 

STANDARD 

DEVIATION 

1.58 

1.61 

COEFFICIENT 
OF  VARIATION 

0.468 

0.485 

TIME  OF 
MAXIMUM 
SIGNAL 

IN  THE 
ENVELOPE 
(msec) 

MEAN 

31.87 

31.87 

STANDARD 

DEVIATION 

0.34 

0.35 

COEFFICIENT 
OF  VARIATION 

0.011 

0.011 

Table  1.  A  comparison  of  some  first  order  statistics  between 
the  real  4fo  and  synthetic  4f0  sampling  frequencies. 
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Figure  A3.  Enlargement  of  time  dependent  phase  difference  case  from  Figure  Al-d. 
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Figure  A4.  Enlargement  of  the  no  time  dependent  phase  difference  case  in  Figure  A2-d. 
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